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FOREWORD 


In  the  determination  of  the  spherical  harmonics  of  the  earth's  gravity 
field  on  the  basis  of  observations  of  artificial  earth  satellites,  the 
separation  of  the  various  terms  in  the  series  is  highly  dependent  on  the 
distribution  of  the  orbital  inclinations  of  the  satellites  used  in  the 
solution.  On  occasion,  one  wishes  to  know  at  what  orbital  inclinations 
it  would  be  most  useful  to  launch  an  additional  satellite  to  improve  the 
gravity  solution  based  on  available  data.  A  complete  computer  simulation 
to  conduct  such  a  study  is  prohibitively  expensive.  This  report  describes 
a  computer  program  to  approximate  the  simulation  and  presents  an  example 
of  such  a  study. 
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INTRODUCTION 

The  coefficients  of  the  expansion  of  the  earth's  gravity  field  of  decree, 
and  order,  m,  are  frequently  divided  into  three  classes 

(a)  zonal  coefficients,  for  which  m  =  0 

(b)  resonance  coefficients,  for  which 

( 1 )  m  a.  k  x  period  of  earth’s  rotation 

orbit  period 

where  k  =  1 ,  2 

(c)  the  remaining  coefficients 

Zonal  and  resonance  coefficients  constitute  a  relatively  small  proportion  of 
the  total  set  of  coefficients;  and,  since  they  nave  very  large  effects  on  the 
motion  of  artificial  satellites,  they  are  relatively  easy  to  determine  to  the 
accuracy  desired  for  non-satellite  applications  such  as  the  computation  of  a 
gravity  anomaly  or  geoid  height  on  the  earth's  satellite.  Therefore,  in  this 
study  attention  was  focused  on  "the  remaining  coefficients",  although  resonance 
coefficients  were  included  as  a  special  case  in  the  computer  program  designed 
to  study  non-zonal  gravity  coefficients. 

The  orbital  perturbations  due  to  gravity  coefficients  are  given  by  Kaula. 
To  zeroth  order  in  orbital  eccentricity  the  perturbation  in  mean  anomaly  is 
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where  <j>  =  (Z-2p)u  +  (Z- 2p)  M  +  m(fi-9) 
S£mpq  =  the  integral  of  Stapq 
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p  =  central  gravitational  constant 
ag=  earth's  semi-major  axis 
n  =  satellite  mean  motion 
a  =  satellite  orbit  semi-major  axis 
io  =  argument  of  per^ee 
M  =  mean  anomaly 
fi  =  right  ascension 
9  =  longitude  of  oreenwich 
t  1  p 

k  =  integer  [(Jl-m)/p] 


c  =  all  values  for  which  binomial  coefficients  (  )  are  non-zero 
q  =  0  for  low  eccentricity 
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The  formulation  and  study  discussed  in  this  report  is  based  totally  on  the 
evaluation  of  this  perturbation  for  reasons  given  below. 


FORMULATION 

The  effect  on  mean  anomaly,  AM,  is  a  maximum  when  the  divisor, 

d  *  ’  H 

i  =  (Cl— 0 )  [(£-2p)  +  (£-2p)  +  m]  ,  is  a  minimum.  The  first  term  is 

small  compared  to  the  other  two  since  the  ratio  of  perigee  rate  to  the  motion  of 

longitude  of  the  node,  j/('s  -  0),  is  small.  The  divisor  is  minimized  either 

for  p  =  1/2  for  any  m  or  for  p  =  £+1  when  equation  (1)  is  satisfied. 

2 

Consideration  of  $  shows  that  the  first  condition  gives  rise  to  the  "m-dailv" 
effects  while  the  second  gives  rise  to  the  resonance  effects:  For  p  =  ?/2  and 
even  degree  coefficients,  the  effect  is  multiplied  by  sines  and  cosines  of 
if  -  m(fJ  -  0),  (neglecting  the  term  in  m)  so  that  the  effects  have  m-daily  periods 
of  (the  earth's  rotational  period/m);  for  p  =  9./ 2  and  odd  degree  coefficients 
=  M  +  m  (f!  -  0)  ,  so  that  the  effect  is  fundamentally  at  the  orbit  period  but 
modulated  at  the  m-daily  period.  For  p  =  £+1,  and  m  odd,  <j>  =  mw  -  M,  which  is 

~~r~ 

the  phase  of  the  resonance  frequency  between  mm  and  M;  for  m  even,  A  =  >!  +  (mm  -  M) 
so  chat  the  resonance  frequency  modulates  the  effect  which  is  fundamentally  at 
the  orbit  frequency. 

The  computer  program  developed  to  study  the  effect  of  satellite  inclination 

on  non-zonal  gravity  parameter  solutions  is  based  on  the  assumption  that  such 

solutions  are  primarily  dependent  on  the  determination  of  the  phase  and  amplitude 

of  these  m-daily  and  resonance  effects  on  mean  anomaly.  For  strict  m-daily 

(even  degree)  and  strict  resonance  (odd  degree)  terms,  this  assumption  is 
\ 

reasonable  at  low  orders  because  the  effects  on  mean  anomaly  are  an  order  of 
magnitude  larger  than  those  on  other  orbital  elements  as  a  result  of  the  squaring 
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of  the  small  divisor.  For  tiie  effects  modulated  at  the  m-daily  (odd  degree) 
period  or  resonance  (even  degree)  period  and  for  higher  orders  or  higher 
eccentricities,  the  effects  on  other  orbital  elements  are  comparable  to  those  on 
mean  anomaly.  For  these  cases,  the  assumption  is  conservative  but  should 
provide  a  guide  as  to  the  relative  usefulness  of  an  additional  satellite  at 
various  inclinations. 

It  is  assumed  that  the  determination  of  the  phase  of  the  effect  serves  to 

separate  the  "C"  and  "S"  coefficients  of  a  given  degree  and  order  so  that  only 

the  amplitude  of  the  effect  need  be  evaluated.  It  is  assumed  that  the  modulated 

and  unmodulated  effects  can  be  separated  so  that  odd  and  even  degrees  are  treated 

separately.  The  problem  then  is  to  form  normal  equations  by  evaluating  the 

partials  of  the  amplitude  of  the  mean  anomaly  with  respect  to  either  odd  or  even 

degree  coefficients  of  a  given  order  which  are  given  by  equation  (2)  with 

S  /  *.  — 

£mpq  dT  ^ 

86400  T  a  10_5  f  . ,  .  .  ]  ‘2 

_ e _  2  ( C-m)  !  (2  P  +  1) 

2tt  - - - —  +  m  L  J 

where  T  is  the  orbit  period  in  minutes,  a  is  the  earth's  radius  in  minutes, 

q  =  i/2  for  diurnal  effects  and  q  =  ({  +  l)/2  for  resonance  effects.  The  10 

_8 

tactor  expresses  the  result  as  10  1  units  in  normalized  coefficients  per  meter 
accuracy  in  the  determination  of  tiie  amplitude  of  the  effect  (that  is,  the 
square  root  of  the  observation  weight  is  included  in  the  partial  derivative). 

For  more  exact  calculation  of  resonance  effects,  the  1436  minutes  should  be 
replaced  by  the  period  of  earth's  rotation  relative  to  the  line  of  nodes  of  the 
satellite;  but  the  approximation  does  not  affect  separation  of  coefficients. 

The  normal  matrices  were  made  for  ten  degrees  at  each  order,  once  for  odd 
degree  coefficients  and  again  for  even  coefficients.  The  equations  are  formed 


for  a  given  set  of  N  orbital  inclinations  plus  an  (N  +  1 )  inclination  which  is 
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varied  over  the  range  of  inclinations  of  interest.  The  eigenvalues  of  the 
normal  matrices  are  computed  and  their  square  roots  are  plotted  versus  the 
inclination  of  the  (N  +  1)  ST  satellite.  Since  the  inverse  of  the  square  root 
of  the  eigenvalue  is  the  standard  deviation  of  a  linear  combination  of  normalized 
gravity  coefficients,  the  optimum  orbital  inclination  for  a  new  satellite  is 
that  inclination  for  which  the  eigenvalues  are  a  maximum. 


EXAMPLES 

Satellites  with  Doppler  transmitters  have  been  launched  at  about  a  dozen 
orbital  inclinations.  However,  many  of  the  satellites  were  launched  many  years 
ago  before  satellite  oscillators  and  ground  equipment  reached  its  present  state 
of  development.  Modern  data  have  been  obtained  at  four  orbital  inclinations, 
63.7°,  90°,  108°,  and  115°.  The  program  described  was  executed  to  determine 
what  satellite  inclination  would  contribute  most  to  the  improvement  in  the 
gravity  coefficients  determined  by  the  set  of  five  satellites.  Since  five 
satellites  would  determine  five  sets  of  even  or  five  sets  of  odd  coefficients 
at  most  under  the  assumptions  used,  only  the  3rd,  4th  and  5th  eigenvalues  of  the 
solutions  were  plotted,  even  though  10  even  order  and  10  odd  degree  coefficients 
were  computed  for  each  order.  Results  for  orders  1-4  and  17  based  on  m-daily 
effects  and  orders  13  and  14  for  resonance  effects  are  given  in  Appendix  A  for 
even  degree  coefficients  and  in  Appendix  B  for  odd  degree  coefficients. 

The  scale  on  the  right-hand  side  gives  the  standard  deviation  of  a  linear 

_g 

combination  of  normalized  coefficients  in  units  of  10  per  meter  accuracy  in 
the  determination  of  the  amplitude  of  the  m-daily  or  resonance  effect.  For  the 
even  degree  low  order  coefficients,  there  is  significant  improvement  in  the 
accuracy  of  the  4th  or  5th  eigenvalue  with  the  addition  of  data  from  a 
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satellite  with  an  inclination  as  little  as  1  or  2  degrees  separated  from  the 
inclinations  and  from  180°  minus  the  inclinations  of  the  first  four  satellites. 
For  the  seventeenth  order  coefficients  and  for  non-resonance  odd  degree 
coefficients,  the  standard  deviations  in  the  coefficients  per  meter  accuracy  of 
determination  of  the  amplitude  of  the  m-daily  effects  are  equal  to  or  larger 
than  the  expected  size  of  the  coefficients.  However,  the  results  are  pessimistic 
due  to  the  neglect  of  perturbations  to  orbital  elements  other  than  mean  anomaly, 
due  to  neglect  of  information  arising  from  orbit  eccentricity,  and  because 
averaging  large  quantities  of  data  will  yield  accuracies  in  amplitude  better 
than  one  meter.  In  view  of  the  pessimism,  the  graphs  lead  to  the  same  conclusion 
as  those  for  even  degree  low  order  coefficients. 

Although  the  charts  show  there  is  a  wide  selection  of  orbital  inclinations 
which  would  improve  the  four  satellite  solutions,  the  determination  of  the 
amplitude  and  phase  of  the  m-daily  effect  can  best  be  done  for  high  inclination 
satellites  because  more  ground  stations  will  observe  higher  inclination 
satellites,  and  because  high  latitude  stations  will  observe  more  passes  per 
day.  Therefore  the  optimum  inclination  for  an  additional  satellite  under  the 
assumptions  of  this  example  would  be  near  either  85°  or  95°. 

Since  96.6°  is  the  inclination  used  by  LANDSAT  type  satellites,  a  further 
experiment  was  run  in  which  it  was  assumed  that  data  were  available  from  five 
satellites  at  inclinations  of  63.7,  90.0,  96.6,  108.0,  and  115.0°.  Appendix  B 
gives  the  accuracy  to  be  expected  in  first  through  fourth  order  coefficients 
as  a  function  of  the  inclination  of  a  sixth  satellite.  Only  two  of  the  graphs 
show  sensitivity  because  for  a  six  satellite  solution,  the  sixth  eigenvalue 
should  have  been  added  to  the  plot.  The  first  and  third  order  even  degree 
results  do  show  sensitivity  because  for  these  coefficients  the  63.7  degree 
inclination  and  the  65°  retrograde  satellite  give  nearly  redundant  information. 
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SUMMARY 


A  computer  program  was  developed  to  help  assess  the  sensitivity  of  solutions 
for  gravity  coefficients  from  satellite  observations  to  the  orbital  inclination 
of  the  satellites  observed.  Sample  cases  indicate  significant  improvement  in 
accuracy  of  the  gravity  field  can  be  achieved  by  the  addition  of  observations 
on  satellites  at  inclinations  separated  by  about  2  degrees. 
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APPENDIX  A 


EIGENVALUES  OF  EVEN  DEGREE  GRAVITY 
COEFFICIENTS  FOR  5  SATELLITE  SOLUTIONS 


A-1 


IGENVALUES  FOR  5  SATELLITE  SOLUTION 
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FIGURE  A-2 .  SECOND  ORDER  EVEN  DECREE  EIGENVALUES  FOR  5  SATELLITE  SOLUTIONS 


GENVALUES  EOR  5  SATELLITE  SOLUTIONS 


FIGURE  A- 3 .  THIRD  ORDER  EVEN  DEGREE  EIGENVAI  UES  FOR  5  SATELLITE  SOLUTIONS 
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FIGURE  A-6 .  FOURTEENTH  ORDER  EVEN  DEGREE  EIGENVALUES  FOR  5  SATELLITE  SOLUTIONS 
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APPENDIX  B 

EIGENVALUES  OF  ODD  DEGREE  GRAVITY 
COEFFICIENTS  FOR  5  SATELLITE  SOLUTIONS 
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EIGENVALUES  EOR  5  SATELL 
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INCLINATION  OF  SATELLITE 
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FIGURE  B-5 .  THIRTEENTH  ORDER  ODD  DEGREE  EIGENVALUES  FOR  5  SATELLITE  SOLUTIONS 


FIGURE  B-6 .  FOURTEENTH  ORDER  ODD  DEGREE  EIGENVALUES  FOR  5  SATELLITE  SOLUTIONS 
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APPENDIX  C 

EIGENVALUES  FOR  6 
SATELLITE  SOLUTIONS 


C-I 


*&&**&*&*<.-.. -uv  ■  - 


j 


EIGENVALUES  FOR  6  SATELLITE  SOLUTIONS 
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FIGURE  C-l.  FIRST  ORDER  EVEN  DECREE  EIGENVALUES  FOR  6  SATELLITE  SOLUTIONS 


EIGENVALUES  FOR  6  SATELLITE  SOLUTION 
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INCLINATION  OP  SATELLITE 


EIGENVALUES  TOR  6  SATELLITE  SOLUTIONS 

ORDER-  4  EVN  DEGREE 


FIGURE  C-4.  FOURTH  ORDER  EVEN  DEGREE  EIGENVALUES  FOR  6  SATELLITE  SOLUTIONS 
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FIGURE  C-5 .  FIRST  ORDER  ODD  DEGREE  EIGENVALUES  FOR  6  SATELLITE  SOLUTIONS 
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Washington,  DC  20350  2 

Naval  Oceanographic  Office 

Bay  St.  Louis,  MS  39522  2 

Office  of  Naval  Research 
Physical  Sciences  Division 
800  N.  Quincy  St. 

Arlington,  VA  22217  2 


U.  S.  Geological  Survey 

ATTN:  Paul  Needham 

Reston,  VA  22091  2 

National  Aeronautics  and  Space 
Administration 
ATTN:  ERG- 2 

Washington,  DC  20546  3 

Jet  Propulsion  Laboratory 
ATTN:  Dr.  William  Melbourne 
Pasadena,  CA  91103  3 

Defense  Mapping  School 

Ft.  Belvoir ,  VA  22060  3 


Copies 

Library  of  Congress 

ATTN:  Gift  and  Exchange  Division 

Washington,  DC  20540  4 

National  Aeronautics  and  Space 
Administrat ion 

Scientific  and  Technical  Library 
Code  NHS  22,  Rm.  BA39 
600  Independence  Ave . ,  SW 
Washington,  DC  20546  2 

Naval  Electronics  Systems  Command 
Navy  Space  Project,  PME106 
Washington,  DC  20360  3 

Office  of  Chief  of  Naval  Operations 
Naval  Oceanography  Division  (NOP-952) 
Bldg.  1,  U.  S.  Naval  Observatory 

Washington,  DC  20390  2 


Naval  Research  Laboratory 
ATTN:  Mr.  A1  Bartholomew  ] 

Dr.  Peter  Vogt  1 

Washington,  DC  20375  2 

Air  Force  Geophysics  Laboratory 
Hanscom  Field 

Bedford,  MA  01731  2 


National  Oceanic  and  Atmospheric 


Administration 
National  Ocean  Survey 
ATTN:  Dr.  John  Bossier  1 

Dr.  William  Strange  1 

Mr.  Larry  Hothem  1 

Rockville,  MD  20850  2 

Goddard  Space  Flight  Center 
ATTN:  Dr.  David  Smith  1 

Mr.  James  Marsh  1 

Dr.  Joseph  Siry  1 

Greenbelt,  MD  20771  2 
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\raval  Space  Surveillance  System 


Dahlgren,  VA  22448  1 

U.  S.  Navy  Astronautics  Group 
Point  Mugu,  CA  93041  2 

Naval  Postgraduate  School 
ATTN:  Prof.  Chris  Mooers  1 

Prof.  Joseph  VonSchwind  1 

CDR  D.  E.  Puccini  1 

Monterey,  CA  93940 


Headquarters  Space  Division  (AFSC) 
Los  Angeles  Air  Force  Station 
Box  92960  Worldway  Postal  Center 


Los  Angeles,  CA  90009  1 

Applied  Research  Laboratory 
University  of  Texas 
ATTN:  Dr.  Arnold  Tucker  1 

Dr.  George  Born  1 

Dr.  Bob  Schutz  1 

Austin,  TX  78712  4 

Department  of  Astronomy 
University  of  Texas 
ATTN:  P.  J.  Shelus  1 

Austin,  TX  78712 

Applied  Physics  Laboratory 
Johns  Hopkins  University 
ATTN:  Harold  Black  1 

Vincent  Pisacane  1 

Johns  Hopkins  Road 
Laurel,  MD  20810  1 


Copies 

U.  S.  Naval  Observatory 
ATTN:  K.  L.  Klepczynski  1 

D.  D.  McCarthy  1 

Washington,  DC  20360 

Pacific  Missile  Test  Center 
Point  Mugu,  CA  9304!  1 

U.  S.  Army  Engineers 

Topographic  Laboratory 

Ft.  Belvoir ,  VA  22060  1 

DMA-IACS  Cartographic  School 
P.  0.  Drawer  934,  Fort  Clayton 
APO  Miami,  FL  34004  5 

National  Geodetic  Survey 
Gravity,  Astronomy  and  Satellite  Branch 
ATTN:  William  E.  Carter  1 

Rockville,  MD  20052 

Institute  for  Laboratory  Astrophysics 
University  of  Colorado 
ATTN:  Dr.  Peter  Bender  1 

Boulder,  CO  80309 

Massachusetts  Institute  of  Technology 
ATTN:  Peter  Morgan  1 

Cambridge,  MA  02139 

Physical  Sciences  Laboratory 
New  Mexico  State  University 
ATTN:  Dan  Martin  3 

Box  3  -  PSL 
Las  Cruces,  NM  88003 


Dept,  of  Earth  &  Space  Science 
University  of  California 
ATTN:  Prof.  William  M.  Kaula  1 

Los  Angeles,  CA  93940 

The  Ohio  State  University 
Dept,  of  Geodetic  Sciences 
ATTN:  Dr.  Richard  Rapp  1 

Dr.  Ivan  Mueller  1 

Dr.  Urho  Uotila  1 

1958  Neil  Ave. 

Columbus,  OH  43210 


The  University  of  Texas  at  Austin 
ATTN:  Dr.  Byron  Tapley  3 

Austin,  TX  78712 
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